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NOMENCLATURE 


C * seal clearance along centerline 

F - axial force 

— 2 2 
F - nondimensional force, F/ 6 uu)(r Q /C) r Q 

h - film thickness 

- integrals defined in Eqs. (21) and (22) 

I ^ • 1 4 " integrals defined in Eqs. (36) and (37) 

- given by Eq. (26) 

J 2 - given by Eq. (41) 

M x - restoring moment 

— 2 3 

M x - nondimensional moment, M x / 6 m<j( r Q /C) r Q 

I 1 2 - transverse moment 

_ 2 3 

M, - nondimensional moment, M,/ 6 uw( r /C) r 
z zoo 

p - pressure 

R - nondimensional radius, r/r Q 

r - radial coordinate 

r' - radial coordinate of extreme pressure 

x,z - orthogonal axes, see Figure 1 



e - tilt parameters, y^/C 

e - angular coordinate 


u - viscosity 


a) - rotational angular velocity 


i 1 i 



Subs grip te : 


ap - approximate 

1 - at inner radius 

m - at midradius 

- at outer radius 


o 



INTRODUCTION 


Seal leakage and short seal life are common problems found In a host of 
Industrial equipment and In other applications. For this reason seal research 
has been carried out since the early sixties in an attempt to understand the 
mechanism of seal operation. Denny (1) has shown experimentally that hydro- 
dynamic effects in a misaligned radial face seal are the cause for axial for- 
ces and pressures in excess of those theoretically predicted for aligned flat 
seal faces. Following Denny's findings, many Investigators have treated this 
problem and there have been several hypotheses put forth to explain the mecha- 
nisms responsible for the development of the lubricating film pressure that 
acts to eparate the primary seal faces. These hypotheses include surface 
angular misalignment [2], surface waviness [3,4,5], microasperities [6], vapo- 
rization of the fluid film [7], and thermal deformation [8,9]. Examination 
reveals that these theories are of limited use from an operation prediction 
standpoint, and that face seal lubrication theory is still very primitive. 

It is also noted that seal dynamics, which is thought to be of major impor- 
tance, is poorly understood. 

In some recent papers [10,11] an attempt is made to solve analytically 
the dynamic behavior of a radial face seal. In these papers the dynamic res- 
ponse of an angular misaligned seal is considered due to a restoring moment 
that coincides with the angular misalignment vector. This restoring moment, 
which is an important factor in seal stability, is by no means the only mo- 
ment acting on the seal faces, a fact that somehow was overlooked in the 
past. A transverse moment that leads the angular misalignment vector by 90 
degrees and is generated by hydrodynamic effects is pointed out in [12], This 
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transverse moment may be the origin of dynamic Instability and has to be con- 
sidered in any dynamic analysis of a realistic seal model. 

In order to establish a better understanding of radial face seal mecha- 
nism of operation, it is Important to obtain the complete system of forces 
and moments. As a first step the hydrostatic effects in a misaligned seal 
were analyzed [12] and both the axial force and tilting moment were found. 

In this paper the hydrodynamic components of the forces and moments will be 
treated. 

In previous works analytical results were limited to very small angular 
tilts 12], or to an approximated film thickness geometry where radial varia- 
tions were neglected (11]. Also the transverse moment, mentioned before, was 
overlooked. 

It is the objective of this paper to present an analytical solution for 
the hydrodynamic effects in a realistic misaligned seal qeometry. This solu- 
tion covers the complete range of seal misal ign.neit from parallel faces to 
surface touch down. 
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ANAl YSIS 


The Reynolds equation for a narrow seal and Incompressible fluid Is: 


3 

3r 


rh' 


!R) 

3r 


6pwi 


^3h 

30 


( 1 ) 


where the film thickness, h, for a misaligned seal (Figure 1) Is given by: 


h * C ♦ yrcose (2) 

It is shown In (12] that while curvature effects may be neglected, thus allow 
Ing the replacement of r In (1) by the mean radius r m , the film thickness in 
(2) should remain a function of both r and e. This is due to the fact that 
the pressure is strongly affected by radial changes in the film thickness 
along the narrow width of the sealing gap. 

In contrast to the film thickness, its circumferential gradient is al- 
most unaffected by the radius and hence may be approximated by 


it! 

36 


y r s i n0 
m 


Eq. (1) then becomes 

f^h 3 f£) ■ - 6 uu ,r m s(ne 


(3) 


Integrating once we have 


J£ 

3r 


Y r m sine(r-r ') 


( 4 ) 
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where r' is a constant of Integration corresponding to the radius where the 
pressure has an extremum. Integrating once again gives 


| 1 P I 

p » - 6u«yr sine 1-* *-(- r ♦ \ ] ♦ C, (5) 

Y^cojre n 2h* 2 Y h^cosu 1 

with Cj as another constant of integration. The boundary conditions for the 
hydrodynamic case are p»0 at r*r Q and at r«r^. Hence 


or 


±) ♦ r ' 

h ° ^ 
o 


h h. 
o 1 


r' * — ! — 12 — 2-J— - C) 
ycose' h~+TT7 u 


( 6 ) 


Since 


h o h 1 


C 2 ♦ C-r{i" 0 +r\ Jcose + > 2 r tvcos^e 


and 


C(h 0 +h^ = 2C 2 + CY(r () +r i )cose 


we get from (6) 

C(r 0 +r j) + 2Yr Q r.cos0 
r 2C + Y ( r 0 +r -j ) cose 
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or In dimensionless form 


R m ♦ cR.cose 

C i . J2 i 

K 1 ♦ eft cose 

III 


( 7 ) 


where R ■ r/r Q , c» Yr Q /C, ano R m ■ (l‘R.)/2. 

Eq. (7) gives the radius where the radial pressure profile reaches Its 

maximum or minimum. This radius clearly differs from the mean radius R at 

m 

which the maximum pressure is obtained when radial variation In h are neglec- 
ted nn. 

The constant In Eq. (5) is found by equating p to zero at r“r Q . Thus 
the pressure distribution Is 


P 


6uwr 

HI 

Y 


Mil II » || 

[4hc ♦ yr'cose) - hh„] 

cos e h^hp 


(8) 


From (6) we have 


C + yr'cose * 2 


Vi 

v h i 


Hence 


6uu)r, 


m 


sine 

1 


cos e h h 


h-h h +h 

"T^OJTVi ‘ hh o ) 


i 


or 


P 


3p<uCR m esi ne 


(r o -r)(r-r.) 



(91 
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where h m » (h^h^/2. 

In (11) where radial variations In h are neglected and the film thick 
ness expression Is: 


h ■ * C ♦ yr^cose 

m ■ m 


the pressure Is given by 


3uwCRes1n0 , 0 . 

p — 3 — (KV • 


m 


which can be rearranged In the form 


p « 3uwCR csine 
ap m 


(VrKr-r,) 


m 


( 10 ) 


The symbol p a _ is used In (10) to Indicate that the pressure In (11) is only 
ap 

an approximation resulting from the omission of radial variations in h. The 
accurate pressure given by (9) is related to the approximate pressure of Eq. 
(10) by 


n A 2 

p ap' h 


(ID 


The film thickness, h, given by Eq. (2) can be written in the form 


h * C(1 + eRcose) 


( 12 ) 


The mean thickness. h„ is 
' m< 
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h » ■ c <' * ‘"m™") 


(13) 


Hence, from (11) It Is clear that for co$e>0, p dp underestimates the pressure 
at any R>R m and overestimates the pressure at any R^. For cose<0, p ap Is 
an overestimation when R>R m and an underestimation of the accurate pressure 
when R<R m - The approximate pressure p ap Is antisymmetric about the line BB 
which connects the highest and lowest points of the seal (see Figure 1), but 
the ratio h m /h Is symmetric about that line. Hence, the accurate pressure p 
Is from (11) also antisymmetric about line BB. 

In the absence of an hydro*'* 1c component this will mean negative pre- 
ssures In the diverging clearance where dh/de>0. For the narrow seal approxi- 
mation cavitation is assumed to occur In this section of the seal (11) and 
the pressure Is assumed to he zero (half Sommerfeld condition). When the 
hydrostatic pressure component does exist the extent of cavitation is depen- 
dent on tne pressure differential across the seal boundaries. The cavitation 
zone decreases as the sealed pressure increases until a full fluid film con- 
dition Is reached. 

The two extreme cases, namely, the half Sommerfeld condition for cavita- 
ting flow and the full fluid film condition for high pressure seals, will now 
be considered. 


CavitatiriLi Flaw : 

Axial Force : 

In the case of a very small pressure differential across the seal boun- 
daries the axial force is 
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J7‘ 


F ■ I I w prdrde 
'1 


(14) 


Substituting Eqs. (12) and (13) Into Eq. (9) we have 


r n ? 

p - auw^riyi-RKR-Rj) 


cslne 


(1 ♦ eR^coseMl ♦ eRcose)' 


(15) 


Neglecting curvature effects and substituting the pressure given by (15), Eq. 
(14) becomes 


F ■ 3uw(-jjr-) r, 


r p\2 2 ( 1 -R) (R-R^ )es1ne 

m o' r| (1*cRcose) z (UtR m cose) 


dRdt) 


(16) 


The Integration over R can be performed by parts noting that 


^ r (l-R)(R-R i ) - 2(R m -R) 


(17) 


Hence 


rl (1-R)(R-R i ) _ (1*R)(R-R 1 ) 

rJ (HcRcose ) 2 * R ' ' «:co*ell*cRco*ej 


‘1 


9 r 1 R— - R 

+ 2 [ m _ 

ccose D j 1 + eRcose 

K i 


(18) 


The first term on the right hand side of (18) vanishs on both limits of the 
integration. The second term yields 


? r 1 R m- R 

^ I rc «p 

ecose D J 1+eRcose 
R 1 


tn( 1+eRcose) 


ecose ecose 


1 


(ecose) 


tjt[ 1 + eRcose - £n(l+eRcose)J } 
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Substituting Into (16) yields 


F • 


6l,ta,( f )2r m I 

(V 


tsine lncose 
o> (.cose) 3 ' '“V 086 


()-«,) 


t sinn 


| • 

(ecose) £ (l ♦ tR cose) 


)de 


( 20 ) 


Defining the Integrals 


Ii(R.e) * I — ■ - s - 1 n W tn(l + cRcose)de 
J (c cose) J 


( 21 ) 


and 


I 2 (e) 


cslne 

p 

(ccose) (1 ♦ cR m cose) 


de 


and using the substitution u»ccose, we find [13]: 


1 1 ( R ,e ) 


cn( UeRcose) , R R 2 . 1 + eRcose 

2(eco*e)* ?ccose ’ 2 cCOse 


( 22 ) 


(23) 


I 2 (0) 


1 + cR^cose 
ccose 


- R «.n 
m 


1 + eRj^cose 
ccose 


(24) 


Substituting Eqs. (23) and (24) into Eq. (20) we have for the axial force 


F 




( 0 ) 


o 


(25) 
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where 


J,(e) - IjO.e) - I 1 (R 1 .•) - (1-R 1 )I 2 (e) 
and F Is a dimensionless force defined by 


(26) 


6 uw(-^) 2 rp 

1 1 ( 1 • 0) and Are obtained from Eq. (23) by substituting R»1 and R*R , 

respectively. After some algebra J^(e) Is obtained as follows: 

, _ 1 1+ecose 1-R 1 1 1+ecose 

1 2(ccos0) ? Weft'. cosb ’ 57cose * 1 £n T+Tf^cose 

R 2 1+eR.cose 

+ T £n 1+eft cose (27) 

m 


It is noted that J^e) as given In (27) is bounded at e«*/2 and hence, F is 
inte^rable over the Interval 0<e<w. This can be readily sho* by expanding 
tne first logarithmic term of (27) in the form 


tn(l+x) * x - ~ + y- + . . . 


(28) 


2 

As e-^/2, cose-K), *Ke orders of (cose) and higher can be neglected, and by 
(28) we have 


lim [ Tj- «n 

e-*i t/ 2 2(ccose) 


1 -R 

1+ecose i_. 

1+eR^cose ” 2tcose J 



Finally, the dimensionless axial force is by (25) and (27) 
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tn 


l + tR, i i. 1+eR 

rar 1 * * ,n 

1 m 


♦ 


HcR 


1/ m 

"r 1 " rar 

in 


in 


1+eR j 

rar: 



(29) 


Eq. (29) gives the axial force over the complete range of misalignment that 
is, from aligned faces at e«0 to touch down at c»l. 

A simpler expression for F can be derived for small tilts by usirg the 
expansion 





(30) 


With this expansion Eq. (29) can be rearranged in the form 


1-R. 3 5 

F - -jle(l-R 1 ) ♦ y(l- R J) ♦ + -’- 1 

t 


+ tc(1 - R m ) + T <1 ’ R m ) + 


+ R ?l c ( R m- R i) + V< R m- R f> + 


(31) 


Neglecting terms of order e and higher and noting that 


'-*m m V R i 


(32) 


we have for the axial force at very small tilts 
F =1- fO- R f ) + fO-R^O+R^lR* 
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which reduces to 


r-jo-v 3 ^ (33) 

The expression for the axial force given by (33) Is the same as the one ob- 
tained In 111] where very small tilts are assumed. Hence, when dealing with 
small c values the omission of radial variations In h Is an acceptable app- 
roximation. 

When the sealing surfaces come Into contact at e*l, the dimensionless 
axial force In (29) becomes 



Using Eq. (32) and noting that 


(1 -t 2 '/Zc Z 

11m in(l-e) * 0 

e-1 


we have at c«l 


p E .i ■ o-V" * V"T?r> R m < 34 > 


Restoring Moment : 

The restoring moment about the x axis (Figure 1) is 


M 


x 


r" 


rr. 



cosedrde 
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Neglecting curvature effects this becomes 


M 


x 



pcosedrdo 


(35) 


The sign convention In (35) assures that a restoring moment will have a posi 
tive value. Using Eqs. (15), (18), and (19), we have 


M 


x 


, / 0x2 3 1 

r m 7 



, cslne 
[ — — — * 

(ecose) 


tn 


1 ♦ccose 
1+eft . cose 


/ 1 p \ tslno 

K i c cose ( 1 +£^ 0056 ) 


ldo 


Defining the integrals 


I 3 (R.e) 


csine 

(ecose) 


j tn(l+cRcose)de 


and 


I 4 («) 


esinedo 


ecose(l+eRcosG) 

m 


(36) 


(37) 


we find 


, ,r ) . «n(l*eRcosa) + R HcReose 
3 ecose ecose 


(38) 


I 4 (e) = ir> 


1+cRj^cose 

ecose 


(39) 


and the nondimens ional restoring moment becomes 
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’(x) 




(40) 


where 


J 2 (e) » I 3 (l.e) - I 3 (R i .e ) - (1-Rj)I 4 (e) 


(41) 


and 


M 


M 


6uu(-£) 2 r 2 0 


After some algebra, J 2 (e) becomes 


1+eR^ose 


1 /n\ . 1 l+ecose . _ 1+ccose D 1 

^2 e ' tcose in 1+cR.cose in Heft m cose " R 1 fcn 1+el^cose 


(42) 


Expanding the first logarithmic term in (42) by the serie of (28), yields 


11m (-tLt ln 3 1 ” R i 

e-u/2 tCOS0 1+cR i cose 1 


Hence, is integrable over the interval 0<e<w and the nondimensional 
restoring moment is 


1 ,2 11* 1 +cR 

M x £n 7^7 * ^ lrv TF ■ lf> T7r ) 
e 1-e m 


1+cR. 3 

e * 771 C ” 


R, 1+ e R m 
+ — ( tn^ — - k— " £r> 


(43) 


m 


For small e, if we use the expansion 
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( 44 ) 


2 3 

tn(l-x) ■ • (x ♦ ~2 * * •..) 


and the one given by (30), Kj x j can be expanded in the form 


4 3 

Mx ■(- -^(l-R 2 ) ♦ ♦ ...) ♦ fled-ig + ^l-R*) ♦ ...) 

2R, 3 , , , 

♦ ~It(R m -R i ) ♦ Y<V R 1> * •••n, («) 


w.ilch, after neglecting high orders of e, becomes for very small tilts 


*x ■ V <<’- R ,) 3 


(45) 


At » s l «.he nondimens ional restoring moment is 


1+R m 7 

■V" 


(46) 


Transverse Moment : 

The moment about the z axis (Figure 1) is 


f* f r 0 ? 

i pr sinedrde 

O'* r J ^ 


which after neglecting curvature effects becomes 


2 , fn f r o 

*z'i j 0 

o' r' 


psinedrde 


( 47 ) 


Again, by Eqs. (15), (18), and (19) 


15 



"• ■ ‘-It?*; J’'rrrTT "*r» 


o J (ccose)' 


- 0-R 1 ) 


csin 2 e 


•7 

(ccose) (1 ♦cR m cose ) 


-)de 


( 48 ) 


Integrating by parts and using Ij and as defined In Eqs. (21) and (22) 
we have 


f 2 

*n(1+cRcose)de * MR.ebine - I I,(R,e)cosede (49) 

'(ccose) • 1 


and 


J 


esin ede 


(ccose) (l+cR |r| cose) 


^(ejsine - 


I 2 (e)cosede 


(50) 


On both boundaries of the integration (e*0 and e**) sine=0. Hence, substitut- 
ing Eqs. (49) and (50) into (48) yields 


M z “ ” 6uu)(-j£) 2 r^ C I -j ( 1 *e ) - I 1 (R 1 ,6) - ( 1-R^ )I 2 (e) Jcosede (51) 


The sum in the brackets under the integral of Eq. (51) was already found, see 
Eqs. (26) and (27). Thus, the transverse moment in its dimensionless form 
can be written as 




1 

J l ‘ c °se 


in 


1+ccose 

1+cR^cose 


(1-R^ 


- ccose(tn 


1+ecose 

l+eR m cose 

m 


2 

+ R‘tn 


1+cR^cose 
1+cR cose^ de 


(52) 


where 
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Z 6 t r °) 2 r 3 

From [14J we find 

f de • .sin-'(tR) 

Q j cose 

Also, integrating by parts, we have 


r 

cosein(HcRcose) 


o J 


cR 


si.^e 

1+t'Rcose 


do 


(53) 


(54) 


and from 115) 


’ T^fo de ■ ' <’ • A*)*! 


(55) 


Substituting Eqs. (53), v 54 ) , and (55) into Eq. (52) yields after some alge- 
bra 


M 2 «{- 7 [sin* 1 (eR i ) - sin" 1 (e) ] ♦ 1 -R i ) ( 1 -c 2 R^)*‘ 

2c 

+ ^[^(l-c^)** - (l-e 2 )*»))R^ (56) 

For small tilts we can use the approximations 

-1 c 3 R 3 

sin (cR) s tR + 2 3 + ••• 


and 
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Thus, for very small c Eq. (56) becomes 



( 57 ) 


Full Fluid Film : 

Under full fluid film condition the hydrodynamic pressure is antisymmet- 
ric about the line connecting the highest and lowest points of the seal ring 
(line BB in Figure 1). This results in a net zero axial force and eliminates 
tne hydrodynamic component of any restoring moment, M^. However, the trans- 
verse moment, M z> becomes twice its value for the half Sommerfeld condition. 
Thus, by (56) we have for full fluid film condition 

M z -l-^(sin" 1 (cR i ) - sin* 1 (c)] ♦ ^(1-R. Ml-e 2 ^)** 

♦ 7(^0 -c 2 ^)** - (l-c 2 ft)nj (58) 

and by (57) for very small tilts 

M, * TfO- R ,) 3 P* (59) 

The hydrostatic pressure ir. a misaligned seal is symmetric about line BB. 
Therefore, in contrast to hydrodynamic effects, it does not produce any trans- 
verse moment Hence, this moment, which is shifted 90 degrees from the mis- 



alignment vector, is entirely due to hydrodynamic effects. 

Although leakage was not treated in the case of a cavitating flow, it 
should be mentioned that under full fluid film conditions the leakage is 
entirely hydrostatic. This Is due to the antisymmetric nature of the hydro 
dynamic pressure which results In zero hydrodynamic leakage. 
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RESULTS AND DISCUSSION 


Values of the nondimensional parameters f, ff z , and the ratio M z / M x 
for the case of cavltating flow (half Sommerfeld condition) are presented In 
Table I and Figures 2 to 4. These cover the whole range of tilt parameters, 
from c»0 to c«l . 

The axial force and the restoring and transverse moments are greatly In- 
fluenced by the radius ratio, i ^/r Q . As /r Q decreases the force and mo- 
ments increase quite rapidly. At a given radius ratio both the axial force 
and the restoring moment increase with c. An increase in the tilt parameter, 
t = vr Q /C, is a result of a decrease in the seal clearance C or an increase 
in the angle of tilt y. Hence, df/dc and dM x /dc are proportional to the 
axial stiffness dF/dC and to the angular stiffness dM x /dv, respecti vely. 

From Table I and Figures 2 and 3 it is seen that decreases in the radius ra- 
tio are accompanied by increases in the axial and angular stiffnesses. Also 
at a given r./r Q the axial and angular stiffness increase with increasing 
tilt parameter e. 

An interesting result is the ratio of transverse to restoring moment 
M z /M x> Whenever t >0.6 M z is smaller than M x , but, as e decreases, the ratio 
M z /M x increases and becomes larger than 1 for any e<0.6. At e*l, the trans- 
verse moment is only about 10 to 20 percent of the restoring moment but be- 
comes about 8 times larger than the restoring moment at e=0.1. These results 
are true for any r./r Q indicating that the center of pressure is not sensi- 
tive to the radius ratio and, what is more important, that there is a strong 
coupling between transverse moment and angular misalignment immediately after 
the misalignment begins. Such coupling can lead to dynamic instability 
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wobbling of the primary seal^slmllar to a whirl In journal bearings. 

In the case of a high pressure seal where the full fluid film condition 
prevails the hydrodynamic components of the axial force and restoring moment 
vanish but the transverse moment becomes twice its value for the cavltatlng 
flow case. The transverse moment vs. tilt parameter for the full film con- 
dition Is presented In Figure 5. 

It Is Interesting to compare the transverse moment due to hydrodynamic 
effects with the hydrostatic tilting moment due to the pressure difference 
across the seal boundaries. These two moments are 90 deqrees apart from each 
other and therefore the problem of dynamic Instability and wobbling in high 
pressure seals is similar to that occuring In low pressure seals with cavita- 
tion . 

As an example, a seal having the following dimensions and operating con- 
ditions was chosen: 


Outer radius, r.., cm . . 

Radius ratio, r Vi» o 

Mean radius, r. cm 

Seal Clearance, C, cm .... 

Shaft speed, n, rpm 

2 

Fluid viscosity , u, N-seo/m . 

2 

Pressure differential, 5p, N/m 


5 

0.9 
4.75 
0.0025 
1000 
1.72 10 ' 
10 6 


In (12) a hydrostatic tilting moment of 430 N-cm is found for that seal 
when e*l. The hydrodynamic transverse moment at the same tilt parameter and 
full fluid film condition is 536 N-cm, and would be even higher for smaller 
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seal clearances C without affecting the hydrostatic tilting moment. Hence, 
the hydrodynamic transverse moment, which was overlooked In the past. Is a 
significant factor In both low and high pressure seals and clearly plays an 
Important role In radial face seal operation. 
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: igure 2 - Nondimensional force as a funct 
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Figure 4 - Ratio of transverse to restoring moment as a 
function of tilt parameter for various radius ratios 
(cavitating flow). 
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